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ABSTRACT 

TreeSPH simulations of galaxy formation in a standard ACDM cosmology, including star formation, 
effects of energetic stellar feedback processes and of a meta-galactic UV field have been performed, 
resulting in a mix of realistic disk, lenticular and elliptical galaxies at redshift z=0. 
The disk galaxies are deficient in angular momentum by only about a factor of two compared to 
observed disk galaxies for simulations with fairly strong star- bursts in early, proto-galactic clouds, 
leading to "blow-away" of the remaining gas in the clouds. In this respect the present scenario 
is hence doing almost as well as the WDM scenarios discussed by Sommer-Larsen & Dolgov. The 
surface density profiles of the stellar disks are approximately exponential and those of the bulges range 
from exponential to r 1 / 4 , as observed. The bulge-to-disk ratios of the disk galaxies are consistent 
with observations and likewise are their integrated B-V colours, which have been calculated using 
stellar population synthesis techniques. Furthermore, the observed /-band Tully-Fisher relation can 
be matched, provided that the stellar mass-to-light ratio of disk galaxies is M/Lj ~ 0.8, similar to 
what was found by Sommer-Larsen & Dolgov from their WDM simulations and in fair agreement with 
several recent observational determinations of M/Li for disk galaxies. 

The ellipticals and lenticulars have approximately r 1 / 4 stellar surface density profiles, are dominated 
by non disk-like kinematics and flattened due to non-isotropic stellar velocity distributions, again 
consistent with observations. 

Hot halo gas is predicted to cool out and be accreted onto the Galactic disk at a rate of 0.5-1 M Q /yr 
at z=0, consistent with upper limits deduced from FUSE observations of OVI. We have analyzed in 
detail the formation history of two disk galaxies with circular speeds comparable to that of the Milky 
Way and find gas accretion rates, and hence bolometric X-ray luminosities of the haloes, 6-7 times 
larger at z ~ 1 than at z=0 for these disk galaxies. More generally, it is found that gas infall rates 
onto these disks are nearly exponentially declining with time, both for the total disk and the "solar 
cylinder'' . This theoretical result hence supports the exponentially declining gas infall approximation 
often used in chemical evolution models. The infall time-scales deduced are ~5-6 Gyr, comparable to 
what is adopted in current chemical evolution models to solve the "G-dwarf problem" . 
The disk of one of the two galaxies forms "inside-out" , the other "outside-in" , but in both cases the 
mean ages of the stars in the outskirts of the disks are >6-8 Gyr, fairly consistent with the findings 
of Ferguson & Johnson for the disk of M31. 

The amount of hot gas in disk galaxy haloes is consistent with observational upper limits. The globular 
cluster M53 and the LMC are "inserted" in the haloes of the two Milky Way like disk galaxies and 
dispersion measures to these objects calculated. The results are consistent with upper limits from 
observed dispersion measures to pulsars in these systems. 

Finally, the results of the simulations indicate that the observed peak in the cosmic star formation 
rate at redshift z ~ 2 can be reproduced. Depending on the star formation and feedback scenarios one 
predicts either a cosmic star formation rate which decreases monotonically with redshift beyond these 
rcdshifts or a second peak at z~6-8 corresponding to the putative population III, and interestingly 
similar to recent estimates of the redshift at which the Universe was reionized. These various scenarios 
should hence be observationally constrainable with upcoming instruments like JWST and ALMA. 
Subject headings: cosmology: theory — dark matter — galaxies: formation — galaxies: evolution — 
galaxies: structure — methods: numerical 



1. INTRODUCTION 

The formation of galactic disks is one of the most 
important unsolved problems in astrophysics today. In 
the currently favored hierarchical clustering framework, 
disks form in the potential wells of dark matter haloes as 
the baryonic material cools and collapses dissipatively. 
It has been shown (Fall & Efstathiou 1980) that disks 
formed in this way can be expected to possess the ob- 
served amount of angular momentum (and therefore the 
observed spatial extent), but only under the condition 



that the infalling gas retains most of its original angular 
momentum — see also Fall (2002). 

Numerical simulations of this collapse scenario in CDM 
cosmologies, however, have so far consistently indicated 
that, when only radiative cooling processes are included, 
the infalling gas loses too much angular momentum (by 
more than an order of magnitude) and the resulting 
disks are accordingly much smaller than required by 
the observations (e.g., Navarro & Benz 1991; Navarro 
& White 1994; Navarro, Frenk, & White 1995; Navarro 
& Steinmetz 1997 using N-body/SPH codes and Bryan 
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2003 using an Adaptive Mesh Refinement code, Bryan 
1999). This discrepancy is known as the angular mo- 
mentum problem of disk galaxy formation. It arises from 
the combination of the following two facts: a) In the 
CDM scenario the magnitude of linear density fluctu- 
ations a(M) =< (SM/M) 2 > 1 / 2 increases steadily with 
decreasing mass scale M leading to the formation of non- 
linear, virialized structures at increasingly early epochs 
with decreasing mass, i.e. the hierarchical "bottom-up" 
scenario, b) Gas cooling is very efficient at early times 
due to gas densities being generally larger at higher red- 
shifts, as well as the rate of inverse Compton cooling 
also increasing very rapidly with redshift (oc (1 + z)*). 
a) and b) together lead to rapid condensation of small, 
dense gas clouds, which subsequently lose energy and (or- 
bital) angular momentum by dynamical friction against 
the surrounding dark matter halo before they eventu- 
ally merge to form the central disk. A mechanism is 
therefore needed that prevents, or at least delays, the 
collapse of small proto-galactic gas clouds and allows the 
gas to preserve a larger fraction of its angular momen- 
tum as it settles into the disk. Weil, Eke, & Efstathiou 
(1998) have shown that if the early radiative cooling is 
suppressed (by whatever means), numerical simulations 
can indeed yield more realistically sized disks — see also 
Eke, Efstathiou & Wright (2000). The physical mech- 
anism by which cooling is suppressed or counteracted, 
however, was not specified. 

Sommer-Larsen et al. (1999, hereafter SLGV99) dis- 
cussed the effects of various stellar reheating mechanisms 
in more detail using numerical TreeSPH simulations of 
disk galaxy formation in the SCDM cosmology. They 
found that reheating of the Universe resulting from more 
or less uniformly distributed star formation does not lead 
to a solution to the angular momentum (AM) problem, 
but that localized star-bursts in proto-galactic gas clouds 
might: if the star-bursts can blow the remaining bulk 
part of the gas out of the small and dense dark matter 
haloes associated with the clouds, then test simulations 
showed that the gas later settles gradually forming an 
extended, high angular momentum disk galaxy in the 
central parts of a large, common dark matter halo. The 
physics of global gas blow-out processes were considered 
in early calculations by Dekel & Silk (1986) and Yoshii 
& Arimoto (1987), indicating that star-bursts might well 
blow out most of the gas in small galaxies with charac- 
teristic circular speed (defined in this paper as the cir- 
cular speed in the disk at 2.2 exponential scalelengths) 
V c < 100 km/s. More recent, detailed simulations by 
Mac Low & Ferrara (1999) suggest that this global blow- 
out scenario may not work so well in disk galaxies, even 
in small ones: the star-bursts typically lead to bipolar 
outflows of very hot gas perpendicular to the disk of the 
small galaxy, only expelling a minor fraction of the disk 
gas. This, however, is due to the particular geometry of 
a disk galaxy (the thinness and flatness of the disk); in 
a more roundish and bulky, proto-galactic dwarf galaxy, 
one would expect the energetic effects of star-bursts on 
the bulk of the dwarf galaxy gas to be much stronger. 

It has proven difficult to realistically implement ther- 
mal energy feedback in cosmological N-body/SPH codes 
incorporating radiative cooling (Thacker & Couchman 
2000 and references therein; Navarro & Steinmetz 2000; 



SLGV99). As first discussed by Katz (1992) this is 
mainly due to problems in resolving a multi-phase in- 
ter stellar medium (ISM) at the scales of individual star 
burst regions in cosmological simulations using SPH. It 
is also related to the smoothing of discontinuities in the 
thermodynamic variables inherent in the SPH method — 
some improvement on this, at least in relation to shocks, 
may be gained by replacing thermal energy by entropy as 
an independent variable in SPH (Springel & Hernquist 
2002). In recent years, however, some progress has been 
made on the implementation of thermal energy feedback 
in SPH simulations (Gerritscn 1997; Mori et al. 1997; 
Thacker & Couchman 2000, 2001; Springel & Hernquist 
2003 — see also Yepes et al. 1997; Ferrara & Tolstoy 
2000; Lia, Portinari & Carraro 2002; Semelin & Combes 
2002; Governato et al. 2002; Abadi et al. 2003 and Meza 
et al. 2003). 

As was shown by Sommer-Larsen & Dolgov (2001, 
SLD01) by going from the CDM structure formation sce- 
nario to the warm dark matter (WDM) one can alleviate 
(and possibly even completely overcome) the AM prob- 
lem without invoking effects of energetic feedback in the 
simulations at all (this does not imply, of course, that 
feedback necessarily is unimportant; it should rather be 
seen as a minimal assumptions approach). Fine-tuning 
of the warm dark matter particle mass to about 1 keV is 
required, however. In contrast, the salient feature about 
"conventional" CDM is that as long as the dark matter 
particles are much heavier than one keV, the actual par- 
ticle mass does not matter for structure formation (note 
though, that axions, despite being ultra- light, behave like 
CDM). 

In this paper we show how a mix of realistic disk, 
lenticular and elliptical galaxies 1 can be obtained in fully 
cosmological (ACDM), gravity/hydro simulations invok- 
ing star formation, energetic stellar feedback processes in 
early, proto-galactic clouds and a meta-galactic UV field. 
This is achieved by treating (albeit in a coarse way) the 
gas in regions where star-bursts are returning energy to 
the ISM as a two-phase medium, consisting of a "cold" 
(T ~ 10 4 K) and a hot (T - 10 6 K) component. Related, 
recent work includes Thacker & Couchman (2001), Stein- 
metz & Navarro (2002) and Springel & Hernquist (2003). 

In section 2 the implementation of stellar feedback pro- 
cesses in the simulations is discussed. Section 3 gives a 
short presentation of the numerical code and the initial 
conditions. The simulations themselves are briefly de- 
scribed in section 4, and the results obtained are ana- 
lyzed and discussed in section 5. Finally, in section 6 we 
summarize our conclusions and present a brief outlook. 

2. STAR FORMATION AND STELLAR FEEDBACK 

The star formation efficiency (SFE, defined as the ratio 
between the dynamical timescale and the star formation 
timescale — see below) in the Galactic disk is quite small 
at present, at most a few percent (e.g., Silk 1997). After 
extensive experimentation we have found that such a low 
SFE in combination with a "conventional" initial mass 
function (IMF) produces too low an energy feedback rate 
in the early and dense proto-galactic clouds to drive the 
remaining gas out of the (cold) dark matter potential 

1 Pictures of some of the galaxies can be seen at 
http://www.tac.dk/~jslarsen/Hubble_Sequence 
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wells in which they have formed. As a consequence, such 
a low SFE at early times does not lead to a solution of the 
AM problem for the CDM + feedback scenario - see sec- 
tion 5. A similar result was found by Abadi et al. (2003), 
but we note that Governato et al. (2002) found for one 
CDM galaxy formation simulation, that the AM prob- 
lem could be solved (to within about a factor of two) 
using low, early SFE + feedback. In any case, for that 
scenario one is still faced with the problem that the stel- 
lar haloes become too massive (Governato et al. 2002) 
and/or the stellar spheroids too dominant and centrally 
concentrated (Abadi et al. 2003; Meza et al. 2003). More- 
over, Governato et al. find for their CDM simulation that 
the number of stellar satellite galaxies is too large com- 
pared to observations of the Milky Way and M31 (the 
"missing satellites problem" ) . 

To solve the AM problem (to within about a factor 
of two; section 5) as well as the other problems men- 
tioned above (section 5; Sommer-Larsen et al. 2003) we 
find that considerably larger energy feedback rates are 
required at early times. This can be achieved, for exam- 
ple, by invoking at early times 1) a considerably larger 
SFE or/and 2) a substantially more top-heavy IMF. In 
this paper, we shall not attempt to give a detailed, the- 
oretical substantiation of these two possibilities, but a 
few, brief plausibility arguments can be given as follows: 

The stars formed in the early, proto-galactic clouds 
will have low heavy clement abundance, as they will be- 
long to the first generations of stars formed in the Uni- 
verse. Hence a physical motivation for case 1) could be 
related to the fact that the strength of stellar winds de- 
creases with decreasing metallicity (Kudritzki et al. 1989; 
Kudritzki & Puis 2000), possibly with a sharp drop be- 
low metallicities of about 1/100 solar (Kudritzki 2002). 
So for increasing metal abundance, the stronger stellar 
winds may regulate star formation so as to decrease the 
SFE, due to the energy and momentum feedback to the 
star forming, molecular clouds. It is also possible that 
the SFE could depend on the thermal history of the gas 
such that gas which has been heated to temperatures 
well above T crit ~10 4 K (below which atomic radiative 
cooling becomes unimportant) and subsequently cooled 
down again has a lower SFE than gas which has never 
been heated above T cr j t . Physically this could be related 
to a possible multi-phase structure of the former type of 
gas, effects of magnetic fields, destruction of molecules 
etc. 

With respect to case 2) some theoretical work on star 
formation indicates that the IMF of the first generations 
of stars was in fact more top-heavy than present day 
IMFs (Abel, Bryan & Norman 2002; Bromm, Coppi & 
Larson 2002; and references in Chiosi 2000). Padoan 
& Nordlund (2002) find that the high-mass slope of the 
IMF becomes steeper at early times (when the magnetic 
fields in molecular clouds were considerably weaker than 
the B ~ [iG at present), but also that the low-mass cut- 
off of the IMF at the same time increases considerably 
such as to effectively make their proposed IMF more top- 
heavy (or more appropriately "bottom-light") at early 
times (A. Nordlund, 2002, private communication). 

In this paper we consider case 1), but we have no reason 
to believe that case 2) would not work equally well. 

We use two distinctly different star formation and feed- 



back modes depending on the thermal history of the 
star forming gas (as chemical evolution is not yet im- 
plemented in the simulations): 

• "Early" star formation mode 

Cool, dense gas, which has always been cooler than 
T crit ~ 10 4 K is assumed to form stars rapidly, i.e. on a 
timcscale comparable to the local, dynamical time td yn - 
We express this star formation timescale as 

_tjyn = 1 J_ 

where the star formation efficiency e e ~ 1 (we find that 
the specific value of e e used is not critical for the outcome 
of the simulations as long as e e ~ 1; we have used e e = 1 
in this work). Such fast star formation is assumed to be 
triggered when the gas density of an SPH particle exceeds 
a certain critical value, chosen to be nn, = 0.3 cm" 3 ; 
we have experimented with other values and found that 
the outcome of the simulations is very robust to changes 
in this parameter. In a star formation event an SPH 
particle is converted fully into a collisionless star particle 
of the same mass, thereby conserving the total number 
of particles in the simulation. The star particle remains 
a star for the rest of the simulation — non-instantaneous 
recycling is being implemented in a forthcoming version 
of the code. The moment of conversion from SPH to star 
particle is, as is customary, determined by a probabilistic 
approach. 

The triggering of an SPH particle for star forma- 
tion may (A) or may not (B) trigger a burst of self- 
propagating star formation (SPSF) in the cold, dense 
gas surrounding it: in scenario "A" not only the SPH 
particle which gets above the critical density threshold, 
but also its neighbouring cold and dense SPH particles 
with densities above nH. c .iow(< ^h,c) are triggered for 
conversion into star particles on their individual, dynam- 
ical timescales. Such self-propagating star formation is 
observed at present in some star- burst galaxies (e.g., in 
expanding super-shells — see Mori et al. 1997). SPSF is 
likely to have been more common at early times (when 
metallicities were lower and the ISM more homogeneous) 
- see McCray & Kafatos (1987). We stress that the 
above implementation of SPSF is just a simple, schematic 
way of modelling complicated (sub-resolution) physical 
processes. 

For scenario "A" we ran 4 series of simulations with 
riH,cjow= 0.05, 0.1, 0.2 and 0.25 cm~ 3 respectively, cor- 
responding to the conversion of approximately 5, 4, 3 and 
2% of the gas mass initially in the simulation into stars 
at z ^ 5 — 6; we dub this the "population III" scenario. 
The 4 different values of riH,o,iow were used to control the 
strength of the star-bursts to check whether the outcome 
of the simulations was sensitive to the choice of this pa- 
rameter. In particular we wanted to check whether con- 
siderable fine-tuning is required to solve the AM problem. 
Early star-bursts converting more than about 5% of the 
initial gas into stars lead to halo star metallicities which 
arc too large compared to observations (see section 5.2 
and also SLGV99), so this sets an upper limit to the 
strength of the early star-bursts. 

In scenario "B" only the initial SPH particle above the 
critical density threshold is triggered for star formation 
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on the dynamical time scale — we dub this the "popu- 
lation II" scenario. 

When a star particle is formed, it is assumed to rep- 
resent a population of stars born at the same time in 
accordance with a Salpeter IMF. It will hence feed en- 
ergy back to the local ISM: during the first 5 Myr only 
stellar winds are assumed to contribute, during the sub- 
sequent about 35 Myr also (and more importantly) type 
II supernova explosions as well (the lightest stars which 
explode as SNII have a mass of about 8 M and a lifetime 
of about 40 Myr). Only stars with masses greater than 
30 Mq are assumed to contribute significant amounts of 
energy to the ISM through winds, and such stars are as- 
sumed to deposit a total of 10 50 ergs per star over their 
lifetime. Stars with masses greater than 8 M Q are as- 
sumed to deposit 10 51 ergs per star to the ISM as they 
explode as SNII at the end of their lifetime. The energy 
from the "star-burst" is deposited in the ISM as thermal 
energy at a constant rate over the lifetime of the burst, 
which is taken to be the above mentioned about 40 Myr 
(Mori et al. 1997). The burst energy is fed back to the, 
at any given time (<40 Myr), 50 nearest SPH particles 
using the smoothing kernel of Monaghan & Lattanzio 
(1985). Part of the thermal energy is subsequently con- 
verted into kinetic energy (by the code) as the resulting 
"super-shell" expands. While a star particle is feeding 
energy back to its neighbouring SPH particles, radiative 
cooling of these is switched off - this is an effective way of 
modelling with SPH a two-phase ISM consisting of a hot 
component (T ~ 10 6 — 10 7 K) and a much cooler com- 
ponent (T - 10 4 K) — see Mori et al. (1997); Gerritsen 
(1997) and Thacker & Couchman (2000, 2001). Scenario 
"A" above typically leads to complete "blow-away" of the 
remaining gas in the proto-galactic "dwarf" galaxy host- 
ing the star-burst, whereas scenario "B" is much more 
gentle, but still leads to a considerable lowering of the 
density of the remaining gas. 

• "Late" star formation mode 

The gas which forms the disks of disk galaxies through 
smooth, dilute and fairly well ordered cooling flows has 
typically been heated to temperatures well above T cr n 
and has subsequently cooled down and settled onto the 
disk to become potentially star forming. Such gas is as- 
sumed to form stars on a much slower timescale 

, tdyn / n 
t*.l = , \A) 

where e/ -Cl. As discussed previously, such a large 
change of SFE between early and late phases of galaxy 
formation is in our experience required to form disks 
with the densities, radial and vertical structures, cold 
gas fractions, sizes and morphologies of observed disk 
galaxies at present. A value of e/=0.0025 was adopted in 
order to obtain (cold) gas fractions in the disk galaxies 
at z — consistent with observational estimates (e.g., 
Sommcr-Larscn 1996: the gas fraction (by mass) in Sa- 
Sc galaxies is -0.10-0.25 — the mean of our 28 SPSF 
disk galaxy runs (see section 4) is 0.17±0.01%, the dis- 
persion 0.07). This adopted value of ei is somewhat low; 
we expect it to increase when non- instantaneous recy- 
cling of gas is properly taken into account in the simula- 
tions. Star formation through this "channel" is assumed 
to take place down to a quite low density threshold of 



0.01 cm -3 . This, however, is still large enough to ensure 
that the gas has cooled down to a temperature T— 10 4 
K, at which the radiative cooling function is effectively 
truncated. For the (typically disk forming) star particles 
formed through this slow mode feedback was completely 
switched off. This is a crude, but effective way of mod- 
elling in our simulations the results of Mac Low & Fer- 
rara (1999) obtained at much higher resolution, that for 
star-bursts in a well-organized, cold gas disk only a small 
fraction of the feedback energy is deposited in surround- 
ing cold gas. 

Finally we note that eq. (2) implies a star formation 
law in the continuous limit (and above the density thresh- 
old) of the form oc pj g 5 , where p cg is the cold (T~ 10 4 
K) gas density. The exponent of 1.5 is both physi- 
cally and empirically (at z~0) well motivated (Kennicutt 
1998). 

3. THE CODE AND THE INITIAL CONDITIONS 
3.1. The code 

We use the gridless Lagrangian A-body and Smoothed 
Particle Hydrodynamics code TreeSPH described in 
SLGV99. Our TreeSPH code is modelled after that of 
Hernquist & Katz (1989). The spline smoothing kernel 
of Monaghan & Lattanzio (1985) is used throughout for 
smoothing of gas dynamical properties as well as soften- 
ing of gravity, as in Hernquist & Katz (1989). Individual 
time-stepping, implemented in a way similar to that of 
Springel, Yoshida & White (2001), is used in the code. 
The system time-step is always the smallest of all indi- 
vidual time-steps. 

We include radiative gas cooling and heating in the 
simulations. The radiative heating corresponds to a red- 
shift dependent, homogeneous and isotropic UV back- 
ground radiation field produced by AGNs and young 
galaxies. This meta-galactic UV field is modelled after 
Haardt & Madau (1996) — the UV field switches on at 
rcdshift z — 6. The radiative cooling function is that of 
a primordial gas, modified by the effects of the UV field 
as discussed by, e.g., Vedel et al. (1994). The code fur- 
thermore incorporates inverse Compton cooling, which 
is also explicitly redshift-dependent. Star formation is 
incorporated by converting SPH particles into star par- 
ticles by the schemes described in the previous section. 
The sum of the numbers of SPH and star particles is kept 
constant. 

The smoothing length of each SPH particle is adjusted 
so as to keep the number of neighbors close to 50. 

Finally, we have in this work used the shear-free 
Balsara (1995) viscosity, rather than the standard 
Monaghan- Gingold (1983) viscosity used previously. 

3.2. The initial conditions 

Our cosmological initial conditions are based on a 
ACDM model with (Q M , ft A )=(0.3,0.7) and a Hubble pa- 
rameter H = 100/ikms -1 Mpc -1 = 65 km s -1 Mpc -1 , 
resulting in a present age of the Universe of 14.5 Gyr. 
Following Eke, Cole, & Frenk (1996) we normalize the 
spectrum to o^{z = 0) = 1.0, where as customary <r| 
is the mass variance within spheres of comoving radius 
8/i _1 Mpc, extrapolated from the linear regime of per- 
turbation growth. 

First, a cosmological A-body simulation (dark mat- 
ter only) with 128 3 particles, a comoving box length of 
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10 h~ Mpc and starting at a redshift of Zi—39 was run 
using the Hydra code (Couchman ct al. 1995). At z=0, 
gravitationally bound dark matter haloes were identified 
by running a friends-of-fricnds algorithm with linking 
length 0.2. 

Second, 12 dark matter haloes (see section 4) selected 
for detailed galaxy formation simulations were resampled 
and rcsimulatcd: For each dark matter halo a galaxy 
forming sub-volumne was determined by tracing all dark 
matter particles, which at z—0 were inside of 4 r 2 oo (^oo 
being the radius inside of which the mean dark matter 
density is 200 times the critical) , back to the initial con- 
ditions at Zj=39. The initial, co-moving linear extent of 
the galaxy forming sub-volume was ^1.5-3 /i _1 Mpc. At 
Zi, SPH particles were added to the galaxy forming sub- 
volume, one SPH per dark matter (DM) particle, with 
a mass of /j, times the original DM particle mass, where 
fb is the baryonic fraction (see below). The masses of 
the DM particles in the sub-volume were reduced ac- 
cordingly to (1 — fb) times the original mass. In the 
surrounding region of the cosmological simulation (with 
co- moving linear extent of about 10 /i _1 Mpc) no SPH 
particles were added and the dark matter particles were 
increasingly coarsely resampled with increasing distance 
from the galaxy forming sub- volume. This, by now cus- 
tomary, procedure has been described by many authors, 
e.g., Navarro & White (1994); Gclato & Sommer-Larsen 
(1999) and Thacker & Couchman (2000). The galaxies 
were selected to be at least 1 Mpc away from galaxy 
groups and 0.5 Mpc away from larger galaxies at z=0. 
They were also selected such as not to undergo major 
merging events (with companion to parent mass ratios 
of 1:3 or more) since z=l. 

We use a baryonic mass fraction fb = 0.10, consis- 
tent with nucleosynthesis constraints (0.015ft. -2 < Of> < 
0.02/i~ 2 ) and with the observationally determined bary- 
onic fractions in galaxy groups and clusters (e.g., Ettori 
& Fabian 1999). The masses of the SPH (and star) par- 
ticles and the high resolution DM particles in the galaxy 
forming sub-volume are 4.0 x 10 6 and 3.6 x 10 7 h~ 1 M Q , 
respectively. The SPH particles are assigned an initial 
thermal energy corresponding to a temperature Tj ~ 
100 K. Gravitational interactions between particles are 
softened according to the prescription of Hcrnquist & 
Katz (1989), with softening lengths of e s ,* = 1.3 ft _1 kpc 
for the gas and star particles and Idm — 2.8 /i _1 kpc for 
the high resolution DM particles. The gravitational soft- 
ening lengths were kept fixed in co-moving coordinates 
until z=2.9, then constant in physical coordinates. 

To check for effects of numerical resolution, one series 
of disk galaxy simulations was run with half these soften- 
ing lengths, and one additional simulation was run with 
one quarter of these softening lengths and eight times 
higher mass resolution. For these latter simulations the 
softening lengths were kept fixed in co-moving coordi- 
nates until z=6.7, then constant in physical coordinates. 

4. THE SIMULATIONS 

We selected 12 dark matter haloes from the cosmo- 
logical N-body simulation for the detailed galaxy forma- 
tion runs. The masses of these haloes span about a fac- 
tor of 10, with 4000-45000 dark matter particles inside 
r2oo at z~0 and circular velocities at r2oo ranging from 
about 100 to 185 km/s. After resampling the galaxy for- 



mation simulations consisted of 30000-150000 SPH+DM 
particles. We started out by running all 12 galaxy sim- 
ulations using the self-propagating star formation pre- 
scription (mode "A" in section 2) with a lower density 
threshold of nH,edow—0.l cm~ 3 (and n^ e =0.3 cm~ 3 , as 
always). Seven of the resulting galaxies had at z=0 dis- 
tinctly disk galaxy like morphologies and kinematics, the 
remaining 5 were lenticular (SO) or elliptical like. Four 
additional series of simulations were subsequently run 
for the 7 disk galaxies: three using again early SPSF 
with nH, e ,iow—0-05, 0.2 and 0.25 cm -3 and one series 
with fast, early, but non self-propagating star formation 
(mode "B" in section 2). Typically, ^30000 system time- 
steps were required to run from 2^=39 to z=0. To test 
for possible numerical effects we ran one additional (and 
more time consuming) series of 7 disk galaxy simulations 
with n_ff, e .zotu=0.2 cm~ 3 , but with gravitational softening 
lengths of 0.65 ft^kpc and 1.4 ft^kpc for the SPH (and 
star) and dark matter particles, respectively — the out- 
come of these simulations compared fairly well with the 
results of the same simulations with "standard" gravi- 
tational softening lengths. We shall briefly discuss this 
in section 5.9. About 50000 system time-steps were re- 
quired to run to z~0 for these simulations. 

In addition we re-simulated one of the disk galaxies 
at an 8 times higher mass resolution (a very CPU in- 
tensive run). This was achieved by extracting the same 
region from cosmological initial conditions with identical 
parameters, but on a 256 3 mesh rather than the 128 3 grid 
used so far. Identical low-wavenumber Fourier modes 
were used for both grids, up to the Nyquist wavenum- 
ber of the coarser mesh. Above that, additional high- 
wavenumber modes were added to the finer mesh (up to 
its Nyquist wavenumber) to account for the extra small 
scale power in this simulation. 

This procedure resulted in SPH (and star) and high 
resolution DM particle masses of 4.9 x 10 5 and 4.4 x 10 6 
h Mq, respectively. Moreover, we used a four times 
higher force resolution for this simulation corresponding 
to gravitational softening lengths of 0.33 and 0.69 /i _1 kpc 
for the SPH (and star) and high resolution DM parti- 
cles, respectively. The simulation was an SPSF run with 
^ff,e.iotu=0.25 cm~ 3 and the result is compared to that of 
the "standard" resolution simulation of the same galaxy 
(also with nH, e ,low=0-25 cm -3 ) in section 5.9. This run 
required ^80000 system time-steps. 

Finally, to illustrate the importance of having a large 
initial SFE we ran the 7 "normal resolution" disk galaxy 
simulations with a constant, low SFE of e=0.02. They 
were of the non SPSF type and run with the standard 
gravitational softening lengths. 

In total we ran 12+4x7=40 "standard" galaxy forma- 
tion simulations, 7+1 additional at higher resolution and 
7 with a constant, low SFE. The computational costs 
amounted to about 5 years worth of single processor CPU 
time on a SGI Origin 2000 computer. 

5. RESULTS 

All results presented correspond to the final state of 
the simulations at redshift z = unless otherwise ex- 
plicitly mentioned. A general presentation of the results 
of the simulations at higher redshift will be given in a 
forthcoming paper. 
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Fig. 1. — Surface density profiles for two Milky Way like disk 
galaxies, denoted SI and S2 in this paper. Radial coordinate is in 
units of the exponential disk scalelength at z=0. Top curves corre- 
spond to SI; stars are marked by solid line, cold gas (T~10 4 K) by 
dotted. We estimate the (stellar) bulge-to-disk ratio of this galaxy 
to B/D=0.17±0.05. Bottom curves correspond to S2 — the curves 
have been displaced 1 dcx downward for clarity (line-types as for 
SI). S2 has B/D=0.27±0.08. 



5.1. Surface density profiles, specific angular momenta 
and structural parameters 

As mentioned in the previous section, 7 of the 12 galax- 
ies simulated have disk galaxy like morphologies and 
kinematics, with the bulk of the stars on approximately 
circular orbits in a disk, most of the rest of the stars in 
an inner, bulge-like component and finally a small frac- 
tion in a round and dynamically insignificant stellar halo 
surrounding the galaxies. The disk galaxies formed in 
our simulations are hence qualitatively quite similar to 
observed disk galaxies like the Milky Way. Of the remain- 
ing 5 galaxies, two have a minor fraction of the stars on 
nearly circular, disk orbits; we classify these as lenticu- 
lars (SOs) and the remaining three galaxies have no disks 
at all; we classify these as ellipticals — see further below. 

The disk galaxies have approximately exponential stel- 
lar disk surface density profiles and exponential to r 1 / 4 
bulge profiles, all in good agreement with observations. 
Two examples of a disk galaxy stellar surface density 
profile arc shown in Figure 1 (for all stars within 2 kpc 
vertical distance from the disk) - the surface density of 
cold gas (T ~ 10 4 K) is also shown. The stellar profiles 
of the lenticular and elliptical galaxies are not exponen- 
tial, as shown in Figure 2, but approximately follow an 
r 1 / 4 law, as shown in Figure 3. 

The bulges of the disk galaxies are generally confined 
to being within a radius tb ~ 1 — 1.5 kpc from the cen- 
ters. Stellar disk scale lengths were determined by fitting 
exponential surface density profiles to the region of the 
disk outside of the bulge, R> rs- 

Bulgc-to-disk ratios were determined by extrapolat- 
ing the exponential disk profiles, obtained as described 
above, to the center of the galaxies. Using these decom- 
positions (which make no assumptions about the bulge 
surface density profiles) the specific angular momenta 



Fig. 2. — Stellar surface density profiles of the three elliptical 
and two lenticular galaxies. 




r'/" [kpc'/"] 



Fig. 3. — Same as in Fig. 2, but versus r 1 / 4 . 



of the stellar disks were estimated taking explicitly into 
account also the region with overlap between disk and 
bulge. As re is comparable to the gravitational soften- 
ing lengths one has to check how affected the structural 
decomposition parameters are by gravity softening effects 
- we return to this at the end of the subsection. 

Characteristic circular speeds V c for the disk galaxies 
were calculated using the approach of SLD01, but as an 
addition taking into account also the dynamical effect of 
the bulges. 

In Figure 4 we show the "normalized" specific angular 
momenta % = j*/V^ of the final disks formed in all 35 
disk galaxy simulations as a function of V c . As argued 
by SLGV99 one expects J* to be almost independent of 
V c on both theoretical and observational grounds. Also 
shown in the figure is the median "observed" value of 
J*, calculated as in SLGV99 and SLD01 for a Hubble 
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Fig. 4. — Normalized stellar specific angular momenta at 2=0, 
J* = j*/V c , of the galaxies formed in the simulations: Disk galax- 
ies formed in the 4 series of runs with early, self-propagating star 
formation (SPSF) arc shown by filled circles. The lenticulars (SO) 
and ellipticals (formed in the SPSF simulations with nH, c ,low=0.10 
cm~ 3 ) arc shown by star symbols and open triangles, respectively. 
Disk galaxies formed in simulations without self-propagating star 
formation are shown by open circles. The solid line shows the me- 
dian value from the observational data on disk galaxies of Byun 
(1992), obtained as explained in SLD01 (but in this paper assum- 
ing _Ho=65 km/s/Mpc), the dotted and dashed lines bracket the 
1-cr and 2-cr intervals around the median. 



parameter of 65 km/s/Mpc, together with the observa- 
tional 1-cr and 2-ct limits. As can be seen from the figure, 
the specific angular momenta of the stellar disks lie only 
about a factor of three below the observed median. This 
is about an order of magnitude better than what is ob- 
tained in similar CDM simulations not invoking stellar 
feedback processes, as discussed by many authors (e.g., 
Navarro, Frenk, & White 1995; SLGV99). 

Also shown in the figure are the normalized specific an- 
gular momenta of the two lenticular and three elliptical 
galaxies. The "effective" V c for these has been calcu- 
lated using the stellar mass versus V c relation for the 
disk galaxies (Figure 14). The specific angular momenta 
of the E/SOs are about an order of magnitude smaller 
than those of the disk galaxies, broadly consistent with 
observations (e.g., Navarro & Steinmetz 1997). 

For the remainder of this paper we shall be concerned 
with the properties of the disk galaxies. We will analyze 
the E/SOs in detail in a forthcoming paper. 

In Figure 5 we show the normalized specific angular 
momenta of the 28 disk galaxies from the SPSF simu- 
lations versus the dimcnsionless spin-parameter A 2 oo of 
their dark matter haloes (A = J^) 1 / 2 /GM 5 / 2 is eval- 
uated at r 2 oo)- Fall & Efstathiou (1980) proposed that 
the angular momentum of a disk galaxy should be pro- 
portional to the spin-parameter of its dark matter halo, 
and indeed a positive correlation is seen in Figure 5, al- 
though the scatter is considerable (the linear correlation 
coefficient is 0.51, so the probability that J* and A200 are 
uncorrelated is less than 3%). Filled triangles, squares, 
pentagons and circles correspond to riH,e,iow — 0.05, 0.10, 
0.20 and 0.25 cm -3 , respectively. It is clear from the 



Fig. 5. — Normalized disk stellar specific angular momenta ver- 
sus the spin parameter A200 (evaluated at r2oo of the dark mat- 
ter halo) for the simulations with self-propagating star formation 
(SPSF). Filled triangles, squares, pentagons and circles correspond 
to nH,e,low = 0.05, 0.10, 0.20 and 0.25 cm -3 , respectively. 




Fig. 6. — As Fig. 5, but for the disk galaxies formed in the 7 
simulations without SPSF (shown by the open circles). Also shown 
(by filled squares) arc the data for the 28 disk galaxies from Fig. 5 
binned into two bins. The vertical bars on these last two data 
points are the dispersions in the two bins; the horizontal show the 
bin size. Finally, shown by star symbols are the results of the 7 
constant, low SFE (e=0.02) disk galaxy runs - no B/D decomposi- 
tion was attempted for these galaxies, because they in general are 
very small and concentrated. 



figure that there is no statistically significant evidence 
that one paxticulcir choice of the para-meter Tifj e i ow leads 
to a better solution of the AM problem than the other. 
Hence, for these SPSF simulations extreme fine-tuning 
of the amount of energetic feedback seems not to be re- 
quired. 

In Figure 6 we show the normalized specific angular 
momenta of the 7 disk galaxies from the simulations 
without SPSF versus A2oo- Also shown are the results 



8 





_ 2 z a _ 
1-ff 

1- CT 

2- a 



200 
[km/s] 




150 200 
V [km/s] 



Fig. 7. — Spin-parameter corrected normalized specific angular 
momenta for the 28 disk galaxies formed in the simulations with 
SPSF. Line symbols as in Fig. 4. 



Fig. 8. — Exponential scale-lengths of all 35 disk galaxies (sym- 
bols as in Fig. 4). Line symbols as in Fig. 4. 



from Fig. 5 binned into two bins. There is clearly an 
indication that the simulations without SPSF arc not 
doing as well in terms of solving AM problem. We also 
show the results of the 7 constant, low SFE (e=0.02) 
disk galaxy runs - these are in general giving very poor 
results in terms of solving the AM problem: the mean, 
spin-parameter corrected (see below), normalized angu- 
lar momentum of these 7 runs is just 27±6% of that of 
the 28 SPSF runs. 

The spin-parameters of our 7 disk galaxy haloes are 
on average somewhat smaller than the median found in 
cosmological N-body simulations — see below. This is 
most likely just an effect of small number statistics. As- 
suming (simplisticly) a linear relationship between J* and 
A 2 oo on the basis of the trend seen in Fig. 5 we correct 
the normalized specific angular momenta of the 28 disk 
galaxies formed in the SPSF simulations as follows to see 
how well we can expect to do in solving the AM prob- 
lem. We assume a theoretical median value of A of 0.05 
(Barnes & Efstathiou 1987; Heavens & Peacock 1988) 
and hence multiply the normalized specific angular mo- 
menta of 28 disk galaxies by (0.05/A2oo)- The results 
are shown in Figure 7. It is seen that the disk galax- 
ies formed in SPSF simulations are really only deficient 
in specific angular momentum by about a factor of two 
relative to the observed median. The fecdback+CDM 
scenario presented in this paper is hence doing almost as 
well as the WDM scenario discussed by SLD01 in terms 
of solving the AM problem for disk galaxies. 

Although we are making an appropriate comparison 
between "observational" disk specific angular momenta 
and model disk specific angular momenta (both obtained 
through B /D decomposition) we note that some authors 
present total (disk+bulge) specific angular momenta for 
simulated galaxies. For comparison with these works we 
calculated j*.tot for the 35 model disk galaxies and found 
< j*.tot >= 0.87 < j* >, so the difference between the 
two types of specific angular momenta is quite small. 

In Figure 8 we show the exponential scale-lengths of 
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Fig. 9. — Spin-parameter corrected exponential scale-lengths of 
the 28 disk galaxies formed in the runs with early, self-propagating 
star formation — see text for details. Line symbols as in Fig. 4. 



the 35 disk galaxies versus V c . Not surprisingly, these 
are also too short by about a factor of two relative to the 
observed median (taken from SLGV99, but corrected to 
ft.=0.65). In Figure 9 we show the same, but corrected for 
spin-parameter effects, as discussed above, and only for 
the SPSF simulations, given the result shown in Figure 6. 

In order to compare our model disk galaxies to obser- 
vations some indicator of Hubble type T is required (as 
usual T=l corresponds to Sa, T—2 to Sab, T=3 to Sb 
and so on). A seemingly obvious choice is the bulge-to- 
disk ratio B / D given the well-known correlation between 
the -B-band bulge-to-disk ratio (B / D) b and galaxy type. 
What we have available from the decomposition of the 
model galaxies are mass B/D. Byun (1992) 2-D decom- 
posed /-band images of ~ 1000 Sb-Sd galaxies (3< T <7) 
from the sample of Mathewson et al. (1992). Contrary to 
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Fig. 10. — Mean /-band bulgc-to-disk ratios versus Hubble Type 
T for the ~ 1000 Sb-Sd galaxies for which Byun (1992) performed 
2-D B/D decomposition. The bars correspond to the dispersion in 
each bin. 



Fig. 11. — Mean birthrate parameter b versus Hubble Type T. 
The original data have been taken from Kennicutt et al. (1994) 
and the data shown have been 3-cr "clipped" . The bars correspond 
to the dispersion in each bin and the dotted line is the linear fit to 
the data used in this paper (eq. 3). 



(B/D)b, /-band B/D trace the mass B/D better (e.g., 
Byun 1992). In Figure 10 we show the mean (B/D)j 
for the ~ 1000 disk galaxies in Byun's sample. The bars 
correspond to the dispersion in (B/D)i for each Hubble 
type. Clearly, it is virtually impossible to classify a disk 
galaxy on the basis of its /-band (and hence approxi- 
mately mass) B/D (at least for 3 < T < 7). 

Another potential disk galaxy type diagnostic is the 
so-called birthrate parameter, b, defined as the ratio 
of the present to past average star formation rate, 
b = SFR/ < SFR >, for the disk stars — see Kennicutt, 
Tamblyn & Congdon (1994) and references therein. In 
Figure 11 we show the mean 6 versus Hubble type (T=l 
7) for the data given in Kennicutt ct al. (1994) after 3-cr 
"clipping". The bars correspond to the dispersion in b 
for each Hubble type. Clearly, & is a much better diag- 
nostic of Hubble type than the mass B/D. A linear fit 
to the data yields 

b = 0.19T-0.13 (1< T< 5) , (3) 

which is the relation we will use in the following when 
comparing our model disk galaxies to observations. 

Figure 12 shows the stellar mass bulge-to-disk ratios 
B/D of the 35 disk galaxies versus the birthrate param- 
eter b. As can be seen from Fig. 11 b ~ 0.1 for Sa's 
increasing to b ~ 1 for Sc's. The lack of trend seen in 
Fig. 12 is consistent with the weak trend found observa- 
tionally for (B/D)j (and also for /T-band B/D ratios, 
e.g., de Jong 1996), when allowing for the large disper- 
sion in each bin seen in Fig. 10. 

As mentioned previously, since tb is comparable to the 
gravitational softening lengths, the structural bulge/disk 
decompositions can be affected by spurious effects of 
gravity softening (e.g., Bate & Burkert 1997; Sommer- 
Larsen, Vedel & Hellstcn 1998a,b). To test for such ef- 
fects we performed the following experiments: Switch- 
ing star formation off and slowly reducing the softening 
lengths to final values half of the initial (which were 1.3 
h~ 1 kpc for the gas (and star) particles and 2.8 h~ 1 kpc 
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Fig. 12. — Bulge-to-disk ratios of all 35 disk galaxies versus 
birthrate parameters — see text for details (symbols as in Fig. 4). 



for DM particles), we continued the 7 disk galaxy SPSF 
runs with Tin e \ ow = 0.20 cm for one additional Gyr, 
starting at z—0. As the orbital time scale in the bulge 
and the main disk is considerably less than 1 Gyr the 
angular momentum distribution of the collisionlcss com- 
ponents (stars and dark matter) will be approximately 
preserved, as angular momenta are adiabatic invariants. 

We then performed the bulge/disk decomposition of 
the resulting 7 disk galaxies as before and found that 
the mean B/D increased by 0.20±0.07 (up from 0.20 to 
0.40), the mean R d by -0.07±0.16 kpc and the mean 
disk j t by 38±73 kpc km/s (up from 229 to 267 kpc 
km/s) as the softening lengths were reduced by a factor 
of two. We do not expect any significant further changes 
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if the softening lengths were reduced even more since now 
£</,* < tb and the central gravitational potential is com- 
pletely dominated by stars and gas. So a statistically 
significant increase in the mean B/D ratio results from 
going to smaller softening lengths, bringing it in better 
agreement with the recent estimate of <0.6 for the mass 
B/D ratio independent of disk galaxy type (Portinari, 
Sommcr-Larsen & Tantalo 2003). As a consequence the 
resulting diskj* would be expected to increase somewhat. 
This is indeed hinted by the above results though the sta- 
tistical significance of this can not be established on the 
basis of the present 7 experiments. The disk scale lengths 
appear to be only slightly overestimated (by 6±14%). 

5.2. Broad band B-V colours of the disk galaxies 

Besides the birthrate parameter, a variety of other 
characteristics of observed galaxies is known to depend 
on their Hubble type (see e.g. the review by Roberts & 
Haynes 1994). In this Section we discuss the integrated 
B-V colours of our simulated disk galaxies. 

We estimate their colours "post-process" : First we cal- 
culate the average chemical enrichment history of each 
galaxy on the basis of its star formation history (SFH), 
inferred from the age distribution of the stars in the final 
galaxy (inside of r=30 kpc) at z = 0. Once the metallic- 
ity evolution Z(t) is computed, we derive from the star 
formation rate ^f (t) of the galaxy its integrated luminos- 
ity L in any given passband as: 



rtG 

-J. 



*(t) L SS p(t G -t,Z{t))dt 



(4) 



where to is the age of the galaxy and Lssp(t, Z) is the 
luminosity of a Single Stellar Population (SSP; Tinsley 
1980; Rcnzini & Buzzoni 1986) of age t and metallic- 
ity Z. In the following we give further details on the 
individual ingredients in calculating the B-V colours. 

5.2.1. The chemical evolution history 

The metal enrichment is calculated post-process, with 
the purpose of obtaining an approximate estimate of the 
integrated B-V colours of the disk galaxies. 

Some level of global pre-enrichmcnt is expected from 
the initial "fast" star formation activity responsible for 
the early feedback effects (see Section 2). This early star 
formation is localized in the cool gas clouds but other- 
wise widely distributed over the simulation volume; it 
typically peaks at redshifts (for the SPSF runs 

- see Section 5.4), followed by a more gentle level of 
activity. A second phase of star formation, associated 
with the formation of the main galaxy itself, commences 
at redshifts z ~ 4 — 5. As a first approximation one can 
hence divide the star formation and chemical enrichment 
history into two phases: an early phase resulting in fairly 
uniform metal (prc)enrichmcnt of the simulation volume 
and a second corresponding to the formation of the main 
galaxy. 

The level of pre-enrichment is calculated with a chem- 
ical evolution model representing a closed-box with the 
same SFH as that of the entire simulation volume. We 
use the chemical evolution model by Portinari, Chiosi & 
Bressan (1998), suitably modified so as to treat the SFH 
*f!(t) as an input information, rather than calculating its 
own star formation rate after a prescribed Schmidt-like 



law. A Salpeter Initial Mass Function (IMF) is assumed 
throughout the calculations, with logarithmic slope 1.35 
and mass limits 0.1 — 100 M & . The metallicity of the 
closed-box at t =1.5 Gyr (z = 4.2) is taken as the level of 
pre-enrichmcnt for the subsequent second phase of main 
galaxy formation. We obtain levels of pre-enrichment in 
the range 1/100 to 1/20 Zq (increasing as the strength of 
the early star-bursts increases, causing a larger fraction 
of gas to be converted into stars during this early phase 
— see Section 2), with 1/40 Z Q as the typical value. 

Starting from the appropriate initial metallicity level, 
we then calculate the chemical evolution of the main 
galaxy based on its SFH during the remaining 13 Gyr 
(from 1.5 to 14.5 Gyr). The SFH is directly provided by 
the age distribution of the star particles that reside in 
the galaxy at the end of the simulation. For the chemi- 
cal evolution of the main galaxy, the closed-box assump- 
tion would be too crude an approximation: not all the 
baryons that end up in the final galaxy are immediately 
and equally available for star formation since the begin- 
ning, as the closed-box model would presume. Realisti- 
cally, what happens is that gas progressively cools out 
and becomes available for star formation (T ~ 10 4 K, 
nu > 0.01 cm -3 ). We can regard this "cool out" as a 
sort of "infall history" for the galaxy. We estimate the 
accretion history of cold gas onto the main galaxy as fol- 
lows. At each time t, we know the total amount of cold 
gas in the simulation volume, M*°*(t), as well as the total 
star formation rate \E' to *(i), and the SFH of the stars that 
end up in the galaxy, ^(t). The rate at which cold, dense 
gas becomes available for the star formation relevant to 
the main galaxy, is simply approximated by: 



M cg (t) = Ml°\t) 



* tot (i) 



(5) 



Hence we calculate the chemical evolution of the galaxy 
assuming an infall model with a mass accretion history 
given by: 

M{t) = M cg (t) + M*(t) (6) 

where M* (t) is the mass of the stars in the final galaxy 

that were formed up to time t. To do so, we further 

adapted the chemical evolution model by Portinari et al. 

(1998) so that it could treat not only the SFH, but also 

the accretion history as input data, in place of the usual 

• _ t_ 
prescription M oc e r . 

Post-process calculation of chemical evolution intro- 
duces some inconsistency with respect to the results 
of dynamical simulations, since the chemical model ac- 
counts for the gas re-ejection from dying stars, which is 
neglected in the simulations. Hence, for the same SFH 
the two models will end up with a somewhat different 
gas versus star content. For the Salpeter IMF adopted 
here, a SSP restitutes ^30% of its mass to the gaseous 
phase over a Hubble time, hence this is the order of mag- 
nitude of the discrepancy between the final mass in stars 
in the chemical model and in the dynamical simulation. 
A consistent calculation of the metallicity distribution 
of the stars in the galaxy would require chemical evolu- 
tion to be implemented in the simulations, so that each 
star particle is consistently labelled with its own age and 
metallicity and can be individually treated as a SSP for 
the sake of colour calculation. This is currently being 



11 



implemented in the code. The above mentioned mis- 
match in final mass in stars, however, mostly affects the 
integrated magnitude of a galaxy, while its colour is only 
marginally sensitive to the slightly different metal enrich- 
ment for different gas fractions. Hence, for the purpose 
of this paper our post-process estimate of the colours is 
adequate. 

5.2.2. The calculation of "integrated" B-V colours 

The chemical evolution calculated above provides us 
with the metal enrichment history Z(t) of the galaxy, 
which combined with its known SFH *f?(t) allows us to 
calculate its integrated luminosity from Eq. (4). As to 
Lssp, we adopt the magnitudes and colours of the SSPs 
by Tantalo et al. (1996). 

The B-V colours, calculated as described above, do 
not include the effects of dust and gas re-processing of 
the stellar radiation, which can be crucial especially for 
the youngest stellar populations. In particular, there is 
evidence from studies of combined dust extinction and 
IR re-emission that the youngest stellar populations can 
be heavily obscured by the parent clouds for periods of 
up to a few 10 7 yr (e.g. Silva et al. 1998; Chariot & Fall 
2000) . In the current paper we deal with this effect in the 
following simple way: We calculate the B-V colour for 
each galaxy by including all stars formed to the present 
(corresponding to no "dust correction" at all) and by 
dropping the light contribution of the stars formed in 
the last 2 x 10 7 yr; we consider the range of colours we 
obtain as our uncertainty and denote the mean of the 
two values (B-V)o- As to the extinction from the dif- 
fuse gas phase, the reddening effect on B-V for spiral 
disks is mild (e.g. Boissier & Prantzos 1999) and more- 
over the observed colours from Roberts & Haynes which 
we use for comparison are the RC3 catalogue "intrinsic" 
(B-V)o colours, which have been corrected for internal 
extinction. 

Further uncertainties and systematic effects on the 
colours may come from the choice of the IMF, which 
influences the level of chemical enrichment. Dynamical 
arguments favour a rather low M/L ratio in spiral galax- 
ies, both for external galaxies and for our Milky Way 
(Sommer-Larsen & Dolgov 2001 and references therein 
- see also section 5.3), implying a more bottom-light 
IMF than the standard Salpeter one adopted here. A 
variety of studies in fact suggest the the slope of the 
IMF gets much shallower than the Salpeter value below 
1 Mq (e.g. Kroupa et al. 1993; Chabrier 2003). In the 
present paper, however, we consistently use a Salpeter 
IMF since our main purpose is to discuss the trend of 
the B-V colours with respect to other galactic proper- 
ties, such as Hubble type etc. 

Finally, systematic differences in the broad-band 
colours of SSPs among different authors exist, while 
better agreement is found for the differential variation 
of SSP colors with age and metallicity (e.g. Chariot, 
Worthey & Bressan 1996; Kodama & Arimoto 1998). 
So, given the above-mentioned sources of uncertainty 
(dust, IMF and theoretical SSPs) the trend of broad- 
band colours with the properties of simulated galaxies is 
to be considered more robust than their absolute values. 

The B-V colours we obtain are compared to the 
observed ones from the review by Roberts & Haynes 
(1994), which is based on a large sample of galaxies com- 
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Fig. 13. — Integrated (B-V)o colours of all 35 disk galaxies, 
calculated using stellar population synthesis techniques — symbols 
as in Fig. 4. Also shown (solid line) are the mean observational 
values for disk galaxies from Roberts & Haynes (1994) and the 
observational 1-a limits (dotted line). 



piled from the RC3 catalogue. We use the disk galaxy 
birthrate parameter b as our morphology indicator and 
convert from the observational morphological type T to 
b as described at the end of Section 5.1 (using eq. 3). 

In Figure 13 we show the predicted (B-V)o colour for 
the 35 disk galaxies versus the b parameter. Also shown 
are the observational results from Roberts & Haynes. 
The agreement between theory and observations is cer- 
tainly satisfactory. 

5.3. The Tully-Fisher relation and mass-to-light ratios 
of disk galaxies 

In Figure 14 we show the stellar mass of the final disk 
galaxies formed in 35 runs versus the characteristic circu- 
lar speed. Also shown is the /-band Tully-Fisher relation 
(TF) of Giovanelli et al. (1997a) for h=0.65, converted 
to mass assuming (stellar) mass-to-light ratios M / Lj = 
0.5, 1.0 and 2.0 in solar units (used throughout) and 
applying an 0.2 mag. offset (Giovanelli et al. 1997b) 
to take into account that the typical disk galaxy in 
our sample is of type Sab. The slope of the "theoreti- 
cal" TF matches that of the observed one very well for 
a constant mass-to-light ratio, which is required to be 
(M/Lj) ~ 0.8. This fairly small value required is consis- 
tent with the findings of SLD01 for their WDM simula- 
tions. Such a low value is fairly consistent with recent dy- 
namically and/or lensing estimated mass-to-light ratios 
for disk galaxies (e.g., Vallejo, Braine & Baudry 2002; 
Trott & Webster 2002), the mass-to-light of the Milky 
Way (SLD01) and can be obtained from stellar popula- 
tion synthesis models provided an IMF somewhat more 
top-heavy (or more appropriately: more bottom-light) 
than the Salpeter law is used (Portinari, Sommer-Larsen 
& Tantalo 2003). The observed chemical enrichment of 
galaxy clusters and the global cosmic star formation and 
enrichment history seem to require an IMF with similar 
properties (Pagel 2002). 
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Fig. 14. — The stellar mass vs. circular velocity "Tully-Fishcr" 
relation for the galaxies from the 35 simulations resulting in disk 
galaxies (symbols as in Fig. 4). Also shown is the observed /-band 
TF relation for Sc galaxies of Giovanclli ct al. (1997) converted to 
mass assuming (M/Li)=0.5 (dashed line), 1.0 (solid line) and 2.0 
(dotted line) and applying an 0.2 mag. offset to correspond to Sab 
galaxies (Giovanclli et al. 1997b), which is the mean type of our 
simulated disk galaxies. Finally, the symbol "MW" with errorbars 
shows the likely range of the stellar mass and characteristic circular 
velocity of the Milky Way (from SLD01, but corrected to stellar 
mass, rather than total baryonic mass). 



5.4. Disk galaxy gas accretion rates 

We have calculated the rates at which hot halo gas 
cools out and subsequently is deposited onto the disks of 
the 35 disk galaxies in our sample. This is of considerable 
interest because observational upper limits can be placed 
on the accretion rate of gas onto the Galactic disk and 
potentially also other disk galaxies. 

As hot halo gas (at T max ~ 10 6 K) starts to cool out 
near the disk and finally is deposited onto it, at some 
much lower temperature (T m i n ~ 10 4 K) highly, but not 
fully, ionized atoms like OVI will be present in the tran- 
sition region. Denoting the accretion column density of 
H atoms per unit time by Nh and the initial (hot halo 
gas) H density nno and neglecting heat conduction, one 
can show that the column density of an ion Z % produced 
between temperatures T m ; n and T max is 



k A z 



< 3 x_ fidT 
( 2 +S) Xe (n H /n H o)A ' [<) 



where k is the Boltzmann constant, Az is the abundance 
of the element relative to hydrogen, s is 1 for isobaric 
and for isochoric cooling, \ lS the number of particles 
per H atom, Xe the number of electrons per H atom, 
fi is the fraction of atoms of element Z which are in 
ionization stage Z % , nu is the H density of the cooling, 
transition-phase gas, and A is the cooling function, such 
that An H n e , where n e is the free electron number den- 
sity, is the energy loss rate per volume (Edgar & Cheva- 
lier 1986, EC). For isobaric cooling (which, cf. EC, is 
the most relevant case for the present problem) EC find 



A(OVI)=3.8-10 14 cm- 2 for N H /n HO =10 7 cm/s and as- 
suming solar gas abundance and composition. 

The result does not change much for hot halo gas oxy- 
gen abundances as low as 1/10 solar: Galactic halo stars 
are observed to be a-element enhanced, reflecting that 
the heavy elements in the stars were primarily produced 
by type II supernovae. For halo stars [O/H]=-1.0 cor- 
responds to [Fe/H]~-1.4 (e.g., Pagel 1997). One would 
expect that the hot halo gas is a-element enhanced as 
well, since the heavy elements in the gas most likely 
originate from a) the early phase of halo star formation 
and/or b) star-bursts in the Galactic disk. In the latter 
case the heavy elements produced by the type II super- 
novae are transported to the halo via fountain flows (e.g., 
Mac Low & Ferrara 1999). For [Fe/H]~-1.4 the radia- 
tive (collisional ionization equilibrium) cooling function 
is almost an order of magnitude less than for [Fe/H]=0.0 
(e.g., Sutherland & Dopita 1993). Hence it follows from 
cq. (7) that A(OVI) should be similar for 1/10 solar 

and solar gas oxygen abundance (for the same Nh/uho)- 
Even if the gas has solar iron-element abundance and is 
a-element enhanced with [O/Fe]~0.4 due to SNII driven 
fountain outflows from the disk the estimate of iV(OVI) 
(at a given Nn/nno) should not increase by more than 
about 50% because of the increased (radiative) cooling 
efficiency resulting from the a-element enhancement. 

From OVI absorption lines in the spectra of ~ 100 
AGNs observed with FUSE (Far Ultraviolet Spectro- 
scopic Explorer), Blair Savage and collaborators find 



A(OVI) sin|6 ~ 1.8-10 1 



and ~ 1.3-10 cm 



m 

the Northern and Southern Galactic hemispheres, respec- 
tively (Savage et al. 2003). Assuming iijjo ~ 10~ 3 cm~ 3 
(EC; B. Savage, 2002, priv. comm.; also what is 
found from our simulations) and combining both sides 
of the disk, this yields a gas accretion rate of 2.9- 10 -3 
M©/yr/kpc 2 locally ("locally" for the observed, OVI ab- 
sorbing gas actually corresponds to a fairly large region 
of the disk out to about 5 kpc from the sun, B. Savage, 
2002, priv. comm.). Assuming a characteristic size of 
the Milky Way's gas disk of i?=15-20 kpc and that the 
above local accretion rate is typical, one obtains an esti- 
mated total accretion rate of about 2.0-3.6 M Q /yr with 
an uncertainty of at least a factor of two. Part (and per- 
haps most) of this is not likely to stem from hot halo gas 
cool-out, but from fountain flows (B. Savage, 2002, priv. 
comm.). 

In Figure 15 we show the present day disk gas accretion 
rates of the 35 disk galaxies versus their characteristic 
circular speed. It is seen that disk galaxies with charac- 
teristic circular speed comparable to that of the Milky 
Way (V c ~ 220 km/s) are found to have accretion rates 
of about 0.5-1 M /yr at z—0, broadly consistent with 
the above observational upper limit. 

As an alternative source of gas infall, in particular in 
relation to the "G-dwarf problem" (see section 5.6), the 
HI high-velocity clouds (HVC) have been considered. We 
note that the accretion rate predicted by our models is 
comparable to the maximum Tosi (1988) could obtain 
with HVCs under a number of favorable assumptions. 
With the "new" source of gas infall discussed in this pa- 
per (cooled-out halo gas) it is no longer required to ap- 
peal to HVCs as the primary source of gas infall onto 
the Galactic disk. 
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Fig. 15. — Disk gas accretion rates at 2=0 for the 35 model disk 
;alaxies (symbols as in Fig. 4). 



Fig. 16. — Normalized disk gas accretion rates per radial bin at 
2=1 and 2=0 for model disk galaxy SI — example of "inside-out" 
disk formation. 



More generally, one can show that the rate at which hot 
gas cools out in the halo is proportional to £x,boi< y >, 
where Lx.boi is the bolomctric X-ray luminosity and 
< A > is the mean, emissivity weighted inverse tem- 
perature of the hot halo gas (Toft et al. 2002). Only few 
relevant X-ray observations of external disk galaxies are 
currently available, but forthcoming XMM-Newton and 
Chandra observations of nearby, nearly edge-on, isolated 
and undisturbed disk galaxies should greatly improve the 
situation and enable the determination of gas accretion 
rates for external disk galaxies from the bolometric X-ray 
luminosity of their hot haloes (Toft et al. 2002). 

For the simulations presented in this paper we have as- 
sumed that the gas has primordial abundance. If the gas 
abundance is 1/3 solar, as in the intra-cluster medium, it 
follows from the work of Toft et al. that the bolometric 
X-ray luminosities and hence gas accretion rates will be 
a factor of 3-4 lower at z=0, so for a Milky Way like disk 
galaxy about 0.1-0.3 M /yr would be expected, even 
more consistent with the upper limit for the Galactic 
disk, deduced from OVI observations. 

5.5. Disk formation: "inside-out" versus "outside-in" 

It is frequently advocated that in hierarchical struc- 
ture formation scenarios, such as the CDM scenario, the 
formation of galactic disks proceeds "inside-out" in the 
sense that an initial formation of the inner parts of the 
disk is followed by a phase of gradual increase of the 
size of the disk lasting for a Hubble time (e.g., White 
& Frenk 1991; Sommer-Larsen 1991; Mo, Mao & White 
1998). With the models presented here we are in a posi- 
tion to test this "inside-out" paradigm. In two of the 12 
dark matter haloes used in this work disk galaxies with 
circular speeds comparable to that of the Milky Way are 
formed. Of the total of eight SPSF simulations run with 
these two haloes, we analyze for the larger halo the one 
with nH,e,iow—0-^5 cm -3 and for the smaller halo the one 
with riH,e,iow—Q-20 cm -3 — the other SPSF simulations 
of these two haloes give similar results, but the runs se- 
lected are the ones which result in the largest j*, being 



19 and 11% larger than the mean for the four SPSF runs 
per halo, respectively. The disk galaxy formed in the 
largest of the two haloes has V c =234 km/s and stellar 
and cold (T ~ 10 4 K) gas masses of M,=5.2-1O 1O M and 
M cg =0.9T0 10 MQ, respectively — we shall denote this 
disk galaxy SI in the following. The other disk galaxy 
formed in the smaller of the two haloes has V c =212 km/s, 
M*=3.5-10 10 M Q and M cg =O.6-lO lo M at z=0 — we de- 
note it S2. These and other physical parameters for the 
two galaxies are given in Tables 1 and 2. 

In Figure 16 we show the rate at which cooled-out gas 
is accreted onto the disk of SI per radial bin at redshifts 
z=l and z=0. The rate has been normalized to unity 
over the entire disk at each redshift and is shown versus 
the radial coordinate in the disk in units of the (stel- 
lar) disk scale length at z—0, Rd- This disk is overall 
clearly forming "inside-out" , but we note that already at 
z=l it is accreting gas out to ~ 6-8Rd, and at z=0 gas 
is still accreting onto the inner parts of the disk. Fig- 
ure 17 is similar to Fig. 16, just for galaxy S2. It is 
apparent from the figure that the formation of this disk 
occurs in a strikingly different, "outside-in" manner. At 
z=l it is, like SI, already accreting gas out to ~ Q-8Rd, 
but at z=Q the gas inflow has contracted to R < 4i?d. 
Hence, though one should obviously not draw any far 
reaching conclusions on the basis of just two Milky Way 
like model disk galaxies, out results indicate that disk 
formation (not surprisingly) is a more complicated pro- 
cess than depicted in the above mentioned analytical or 
semi-analytical works. 

That gas is being deposited onto the disks out to 
~ 6-8i?d is also interesting in relation to the observa- 
tional finding that present day HI disks often extend well 
beyond the (stellar) Holmberg radius i?26.5 (~ 4-5 Rd 
- see van der Kruit & Searle 1982), in fact to about 
1.5i?26.5 (e.g., Wolfe et al. 1986). We shall in a forth- 
coming paper discuss the implications of our models in 
relation to observed medium to high-z damped Lya sys- 
tems (DLAs). 
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Table 1 

Masses 51 of and number of particles' 1 in selected disk galaxies at z=0 



Galaxy 


M* 


M cg b 


M hg c 


Mh,tot d 


N» 


N b 


N hg c 


N b , tot d 




[1O 1O M ] 


[10 10 M Q ] 


[10 10 M Q ] 


[10 10 M Q ] 










SI 


5.16 


0.91 


0.03 


6.10 


8493 


1495 


54 


10042 


S2 


3.46 


0.61 


0.03 


4.10 


5697 


1001 


56 


6754 


S3 


1.48 


0.46 


0.03 


1.97 


2431 


763 


47 


3241 


S3-8 


1.65 


0.20 


0.03 


1.88 


21713 


2684 


381 


24778 



a Within galactocentric distance r=30 kpc 
b Gas with T < 3 ■ 10 4 K 
c Gas with T > 3 ■ 10 4 K 

d Total baryonic mass and number of particles 



Table 2 

Other physical properties of selected disk galaxies at z=0 

Galaxy V c J, jT a ~bJd R~ d M^ s SFR b Type 
[km/s] [kpckm/s] [kpc km/s] [kpc] [M /yr] [M /yr] 

51 234 532±21 1147 0.17±0.05 1.75±0.10 0.54 0.98 0.22 Sab 

52 212 400±23 387 0.27±0.08 1.43±0.07 0.71 1.70 0.62 Sbc 

53 163 236 664 0.28 1.28 0.30 0.22 0.19 Sab 
S3-8 166 248±20 936 0.28±0.13 1.17±0.05 0.16 0.26 0.20 Sab 
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Fig. 17. — Normalized disk gas accretion rates per radial bin at 
z=l and z=0 for model disk galaxy S2 — example of "outside-in" 
disk formation. 



The total gas accretion rates for galaxies SI and S2 (in- 
tegrated over the entire disk) are 0.5 and 0.7 M Q /yr at 
z=0 and 3.0 and 4.5 M Q /yr at 2=1, so the gas accretion 
rates are ~ 6-7 times larger at 2=1 than at 2=0. The 
masses of SI and S2 at z=l are M»=3.5 and 1.8 T0 10 M Q 
and M cg =1.7 and O.9TO 1O M , respectively. Given the 
considerable drop in gas accretion rates from z=l to 2=0, 
that the galaxies are accreting gas out to ~ 6-8Rd at 
z=l and the fairly modest increase in total mass from 



z=l to z=0, it seems safe to conclude that the formation 
of isolated large disk galaxies was well under way by red- 
shift z ~ 1. This seems in line with the observations by 
Simard et al. (1999), which indicate that disk galaxies 
have similar sizes at z ~ 1 and 2=0. 

In Figure 18 we show for galaxies SI and S2 the mean 
age of disk stars (within 2 kpc vertical distance from 
the disk) versus radial coordinate in units of the (stel- 
lar) disk scale length at z—0. The bars correspond to 
the dispersion in each bin. Only stars with formation 
redshift 2/ < 2.6 (corresponding to ages < 12 Gyr) were 
included — this eliminates essentially all halo stars and 
also partly the bulge stars for these two galaxies. In the 
disk of SI there is essentially no age gradient and the 
mean stellar age is ~ 7-8 Gyr corresponding to a mean 
formation redshift of ~ 0.8 (in the two innermost bins 
the mean age is larger due to contribution from bulge 
stars). The disk age dispersion is ~ 2-3 Gyr correspond- 
ing to a dispersion in redshift of ~ 0.3-0.4. At first it 
seems surprising that a disk, forming inside-out, has no 
age gradient. However, firstly the star formation history 
is generally different from the gas accretion history, and 
secondly the star formation rate depends non-linearly on 
the cold gas density, cx pj g 5 (cf. eq. 2). This makes 
the star formation rate quite low in the outer disk at late 
times and hence the average stellar age quite high there 
(but some young stars are present — see below). 

The mean age of the stars in the outer disk of galaxy 
S2 is ~ 6-7 Gyr — slightly less than the outer disk of 
SI. The mean age decreases with decreasing R to about 
4 Gyr in the third innermost bin. This is not surprising 
given the gas accretion patterns at z=l and z—0 shown 
in Fig. 17. In the two innermost bins there is some 
contribution from bulge stars (with Zf < 2.6) increasing 
the mean age there to 6-10 Gyr. 
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Fig. 18. — Mean ages of disk stars in galaxies SI (solid circles) 
and S2 (open triangles). Only stars with formation rcdshift Zf <2.6 
(or equivalently age < 12 Gyr) are included. The bars show the 
age dispersions in each bin. 



Fig. 19. — Distribution of stellar ages in the disks of galaxies 
SI (thick curves) and S2 (thin curves) as a function of the radial 
coordinate in the disk. Only stars with zj <2.6 (or equivalently 
age < 12 Gyr) are considered. The fraction of these younger than 
3, 6 and 9 Gyr are marked by dashed, dotted and solid lines, re- 
spectively. The inner ~ l.bRd are dominated by bulge stars. 



To gain more insight than what can be obtained from 
mean ages, we show in Figure 19 the distributions of 
stellar ages as a function of the radial coordinate in the 
two disks (again for stars with Zf < 2.6). The fractions 
of stars younger than 3, 6 and 9 Gyr are shown for the 
disk of SI with thick curves and for S2 with thin ones. 
For SI the fractions of stars younger than 3 and 6 Gyr 
increase steadily with R out to R ~ 4-5Rd- So the 
distribution of stellar ages does reflect that the disk of SI 
formed inside-out. At R ~ 5-6Rd, the fractions of stars 
younger than 3, 6 and 9 Gyr are ~ 15, 30 and 60%. 

The distribution of stellar ages in the disk of S2 reflects 
the outside-in formation of this disk. At R ~ 5-6i?d, 
there are almost no stars with ages less than 3 Gyr, about 
35% with ages less than 6 Gyr, and only about 10% with 
ages larger than 9 Gyr. The latter explains the slightly 
lower mean age of the outer disk of S2 compared to SI. 

Ferguson & Johnson (2001) find in a recent observa- 
tional study of the outskirts of M31's disk that the typi- 
cal stellar age at ~ 5-6 disk scale lengths is > 7-8 Gyr. 
This seems in fair agreement with our results on the 
mean stellar ages at ~ 5-6 Rd, which are also lower lim- 
its since we have excluded all stars with Zf > 2.6. But at 
least in the disk of SI we find a non- negligible fraction of 
stars younger than 3 Gyr in the outskirts of the disk — 
whether this is consistent with the observational findings 
of Ferguson & Johnson remains to be seen. 

The mean age of the disk stars inside of R = 6Rd, 
with Zf < 2.6 and excluding the two innermost bins, 
which have some bulge contribution, is 7.3 Gyr for SI. 
This corresponds to a mean disk star formation redshift 
of < Zf >> 0.8, so for this galaxy even the formation 
of the bulk of the stellar disk happens at z ~ 1. For 
galaxy S2, the corresponding number is 4.9 Gyr, so for 
this galaxy the stellar disk is arguably still forming — its 



present SFR is also quite respectable, cf. Table 2. 

Finally, our results above on age gradients indicate that 
the colour gradients observed in galactic disks (with the 
disk generally becoming bluer with increasing R) are a 
metallicity effect, rather than an age effect. 

5.6. Gas accretion histories for galactic disks and the 
"G-dwarf problem" 

Observational studies of the metallicity distribution 
of G and K stars in the solar neighbourhood show 
that the distribution is very narrow (peaking around 
[Fe/H]=-0.2 dex; Wyse & Gilmore 1995; Rocha-Pinto & 
Maciel 1996; Hou et al. 1998; J0rgensen 2000; Flynn & 
Morell 1997; Kotoneva et al. 2002), much too narrow 
to be consistent with what is expected from a simple, 
closed-box model of galactic chemical evolution. This 
famous, apparent inconsistency has been coined the "G- 
dwarf problem" and was first discussed by van den Bergh 
(1962), Schmidt (1963) and Pagel & Patchett (1975). 
The currently most popular (and very sensible — see 
below) answer to the G-dwarf problem is to assume that 
the solar neighbourhood was formed gradually over the 
Hubble time through prolonged infall of fairly low metal- 
licity cooled-out gas originating from the hot halo; such 
"infall models" were first introduced by Larson (1972) 
and Lynden-Bell (1975). Assuming, as is often done for 
simplicity, that the infall rate is an exponentially decay- 
ing function of time an infall timescale of ti n f =7— 10 Gyr 
is required to get the models to match recent observations 
(Chiappini et al. 1997; Portinari et al. 1998; Boissier & 
Prantzos 1999). If a gaussian form for the infall rate is 
assumed, then the typical timescale is ~5 Gyr (Prant- 
zos & Silk 1998; Chang et al. 1999). With the current 
models we are able to test these assumptions: 
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Fig. 20. — Rates of infall of cooled-out halo gas onto the disk Fig. 21. — As Fig. 20, but for galaxy S2. 

of galaxy SI versus age. Solid line: all disk, dotted line: "solar 
cylinder" (operationally defined as 2.lR^<R<3.2R l i). 



In Figure 20 and 21 we show the histories of accre- 
tion of cooled-out gas for the disks of galaxies SI and 
S2, respectively. In each figure we show by solid line 
the accretion rate history of the entire disk and by dot- 
ted that of the "solar cylinder" (operationally defined as 
2.1Rd<R<3-2Rd', we assume that the scalelength of the 
Galactic disk is 3.0 kpc, as a compromise between var- 
ious recent determinations, and that the solar circle is 
at i?o=8-0 kpc). For SI the rate of gas infall onto the 
"solar cylinder" is almost perfectly exponentially decay- 
ing with an e-folding time of ~ 5 Gyr (the earliest bin 
in the figures correspond to 2=2.2). For S2 the infall 
history looks more like the one derived by Prantzos & 
Silk (2003) for the solar cylinder using a combination of 
observational data and theoretical arguments. From a 
linear fit to the infall rate history of the "solar cylinder" 
over the age range 0-12 Gyr we derive an infall timescale 
of ~ 6 Gyr. So the infall timescales derived fall short 
of the 7-10 Gyr required by current chemical evolution 
models, but only marginally so. Moreover, we remind 
the reader again that we have been analyzing just two 
Milky Way like disk galaxies in detail. 

5.7. Hot gas in galaxy haloes and pulsar dispersion 
measures 

The amount of hot gas (T ~ 10 6 K) in the Galac- 
tic halo can be constrained from dispersion measures to 
pulsars in the Galactic halo. Of particular interest is 
one pulsar in the high latitude globular cluster M53 and 
three in the LMC (Large Magellanic Cloud) , because the 
distances to these pulsars are known. Rasmussen (2000) 
finds from these data and some modelling of the hot halo 
an upper limit on the cumulative amount of hot gas to 
the distance of the LMC of 1.5-2-10 9 M Q . In Figure 22 
we show the cumulative mass of hot gas vs. galactocen- 
tric distance for three Milky Way like disk galaxies. Solid 
curves are for gas with T > 10 5 K, dotted for T > 10 6 K. 
Top set of curves are for disk galaxy S2, middle for SI 
and lower are for a simulation of SPSF type of the halo 
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Fig. 22. — Cumulative mass of hot gas vs. galactocentric dis- 
tance for three Milky Way like galaxies. Solid curves are for gas 
with T > 10 5 K, dotted for T > 10 6 K. Top set of curves are for 
disk galaxy S2, middle for SI and lower are for a simulation of 
SPSF type of the halo corresponding to galaxy SI, but with a 
gas abundance of 1/3 solar (see Toft et al. 2002), rather than the 
primordial gas abundance used in this work. 



corresponding to galaxy SI, but with a gas abundance of 
1/3 solar (see Toft et al. 2002), rather than the primor- 
dial gas abundance used in this work. It is clear that all 
three hot gas haloes satisfy the above constraint increas- 
ingly so with increasing gas abundance, as discussed by 
Toft et al. 

With the present models we can, however, do better 
than this, because we can "insert" M53 and the LMC 
in our model galaxy haloes and calculate the expected 
dispersion measures directly: 

For electromagnetic waves propagating through a 
plasma, the group velocity v gr will depend on the fre- 
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quency v of the radiation. Specifically, for radio waves 
in an unmagnetized plasma (e.g., Falthammer 1995) 



(1 



u 2 



(8) 



where v pe = (n e e 2 /7rm e ) 1 / 2 is the electron plasma fre- 
quency, e the electron charge, m e the mass of the electron 
and c the speed of light. A pulse travelling a distance L in 
the plasma will have a travel time t given approximately 

by 

L e 2 D m 

* « - + 7i ^ 9 

c 27rm e c^ z 

rather than the "vacuum value" L/c (e.g., Spitzer 1978), 
where the dispersion measure D m is defined as 



D m = 



,dl 



(10) 



Detected radio signals from pulsars show a frequency- 
dependent delay that can be used to evaluate D m . For 
the pulsar in the globular cluster M53, located at a 
distance of «18.9 kpc at (/,&) = (333°,80°) one finds 
D m =24.0 cm -3 pc. For the three pulsars in the LMC, 
located at a distance of ^49.4 kpc at (/, b) ~ (279°,-33°), 
Z3 m =115±13 cm -3 pc is found (Rasmussen 2000). 

To calculate what values of D m we predict to M53 
and the LMC for the haloes of SI and S2 we average 
over the hot halo gas free electron densities in spherical 
shells and correct the densities for resolution effects at 
the disk-halo interface. We calculate D m as the average 
of the value at z=0 and the values obtained at times 100 
and 200 Myr prior to this. One has to adopt a value 
for Ro, the distance of the sun from the galactic center, 
when calculating D m for the model galaxies. Though 
the stellar disks are too small compared to reality, this is 
not necessarily the case for the hot halo. We hence did 
the calculations of D m with two values of i?o=5 and 8 
kpc. For the M53 pulsar we find for SI 24±2 and 16±1 
and for S2 28±1 and 17±1 cm" 3 pc for R =5 and 8 kpc, 
respectively, in fairly good agreement with the observed 
value. For the LMC we find for SI 26±2 and 18±1 and 
for S2 31±1 and 20±1 cm" 3 pc for R Q =5 and 8 kpc, 
respectively. These values are considerably smaller than 
the observed D m — this is probably due to additional 
contributions to the observed D m from the Gum Nebula, 
a largely ionized region at I ~ 230°-290° and distance of 
about 0.5 kpc (Taylor & Cordes 1993), and the diffuse, 
ionized Reynolds-Layer (Reynolds 1993) in the direction 
towards the LMC (Rasmussen 2000). 

It is clear from Fig. 22 and the discussion in Toft 
et al. (2002) that the predicted values of D m above are 
upper limits in the sense that had we included chemical 
evolution in the simulations it would have lead to some 
non-zero level of metal enrichment in the hot halo gas, 
which in turn would have caused a reduction in the pre- 
dicted values of D m . 

5.8. Star formation histories 

As an example of the typical star formation histories 
in our simulations we show (for the same dark mat- 
ter halo) in Figure 23 the star formation rate in units 
of the current one for three of the five simulations re- 
sulting in the largest disk galaxy (with V c — 220-230 
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Fig. 23. — Star formation histories in three of the five of our 
simulations forming the largest disk galaxy (with V c — 220-230 
km/s). Two are for self-propagating star formation simulations 
with n H e ; ot „=0.05 and 0.25 cm -3 , and are shown by dotted and 
dashed lines respectively. The last one is for the simulation without 
self-propagating star formation and is shown by the solid line. 



km/s). Two are with self-propagating star formation 
with riH,e,iow— 0-05 and 0.25 cm -3 and one is for the 
simulation without self-propagating star formation. All 
three reproduce the observed peak in the SFR history 
at 2^2 (Madau et al. 1996), but whereas the simula- 
tions with self-propagating star formation have a low at 
2~4.5-5.5 and another peak at z~6-8 corresponding to 
the putative population III, the one without self prop- 
agating star formation has a monotonically decreasing 
SFR with increasing redshift. Though it would obvi- 
ously be premature to take the star formation histories 
of our simulations of the formation and evolution of in- 
dividual and fairly large galaxies as representative of the 
SFH of the Universe as a whole, our results suggest that 
the global SFH could possibly be bimodal. In relation 
to this it is interesting that recent observations of high- 
z quasars may indicate that the Universe was reionized 
at a redshift of about 6 (Becker et al. 2001; Djorgovski 
et al. 2001; but see also Hu et al. 2002, who advocate 
Zreion ^ 6.6 from observations of a candidate z=6.56 
galaxy; this is, however, still consistent with the location 
of the second SFR peak discussed above). The different 
SFH scenarios should be observationally constrainable 
with upcoming instruments like JWST and ALMA. 

5.9. Effects of numerical resolution 

To check for effects of numerical resolution we did the 
following two tests: 

1. We ran a series of 7 SPSF disk galaxy forma- 
tion simulations with nH, e ,low = ^-2 and with gravitational 
softening lengths equal to half the "standard" values, 
i.e. 0.65 ft. _1 kpc and 1.4 /i _1 kpc for the SPH (and star) 
and active DM particles, respectively. The median nor- 
malized stellar specific angular momentum was 11±15% 
lower than for the "standard" simulations. This was 
mostly due to the median V c being 4±12% larger than for 
the standard simulations. So there may be a small (and 
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statistically insignificant) effect of the gravitational soft- 
ening on the resulting normalized specific angular mo- 
menta - but see below. 

2. The second test was to re-simulate one of our disk 
galaxies, which we shall denote S3 in the following, at 
8 times higher mass and 4 times higher force resolution, 
as detailed at the end of section 4. We shall denote the 
higher mass and force resolution version of S3 by S3-8 
in the following. Masses, numbers of particles and other 
physical characteristics of S3 and S3-8 at z=0 are given 
in Tables 1 and 2. The total (baryonic) masses of S3 
and S3-8 are very similar, differing by less than 5%. The 
stellar mass of S3-8 is 11% larger than that of S3 and 
the cold gas mass correspondingly smaller. We believe 
that it is the combination of the ability to resolve gas 
clouds of smaller masses and the additional small scale 
power that leads to the somewhat higher stellar mass 
and correspondingly lower cold gas mass of the high res- 
olution galaxy S3-8 relative to S3. The characteristic 
circular speed of S3-8 is just 2% larger than that of S3 
despite the 4 times higher force resolution. This may in- 
dicate that force resolution used in our "standard" sim- 
ulations is adequate. The bulge-to-disk ratio of S3-8 is 
B/D=0.28±0.13. To obtain the j* and R d given in Ta- 
ble 2 for S3 we forced the bulge-to-disk decomposition 
to have B/D=0.28 for S3 to enable a fair comparison be- 
tween S3 and S3-8. We find that j« for S3-8 is 5±8% 
larger than that of S3 and that the specific angular mo- 
mentum of the cold gas j cg is 41% larger. The gas accre- 
tion rate is somewhat lower (this is presumably related to 
the larger stellar mass and smaller cold gas mass of S3- 8 
relative to S3) and the star formation rates and birthrate 
parameters, b, are very similar. 

We conclude that the high mass resolution simulation 
generally is in good agreement with the "standard" res- 
olution one, with a tendency for the stellar mass to be 
somewhat larger, the cold gas mass somewhat smaller 
and the specific angular momenta somewhat larger. All 
in all it is hence possible that increased mass resolution 
will be beneficial also with respect to the AM problem. 

6. CONCLUSION AND OUTLOOK 

In conclusion we have obtained a mix of realistic disk, 
lenticular and elliptical galaxies 1 in our ACDM galaxy 
formation simulations with effects of energetic stellar 
feedback processes included. We find that the disk galaxy 
angular momentum problem can be considerably allevi- 
ated (though not entirely solved as yet) in this way, pro- 
vided that the early star-bursts are fairly strong, convert- 
ing 2-5% of the initial gas mass into stars. For early star- 
burst strengths in this range, the improvement on the 
AM problem is about the same, so extreme fine-tuning 
of the feedback seems not required. 

The stellar disks have approximately exponential sur- 
face density profiles and those of the bulges range from 
exponential to r 1 / 4 , as observed. The bulge-to-disk ra- 
tios of the disk galaxies are consistent with observations 
and likewise are their integrated B-V colours, which have 
been calculated in a simplified way ("post-process") us- 
ing stellar population synthesis techniques. Furthermore 
we can match the observed /-band Tully-Fisher relation, 

1 Pictures of some of the galaxies can be seen at 
http: / / www.tac.dk/~jslarscn/Hubble_ Sequence 



provided that the stellar mass-to-light ratio of disk galax- 
ies is M I Li ~ 0.8, similar to what was found by Sommer- 
Larsen & Dolgov (2001) from their WDM simulations 
and in fair agreement with several recent observational 
determinations of M/Lj for disk galaxies. 

The ellipticals and lenticulars have approximately r 1 / 4 
stellar surface density profiles, are dominated by non 
disk-like kinematics and flattened due to non-isotropic, 
stellar velocity distributions, again consistent with ob- 
servations. 

We predict that hot halo gas is cooling out and being 
accreted onto the Galactic disk at a rate of 0.5-1 M /yr 
at z=0, consistent with upper limits deduced from 
FUSE observations of OVI. We have analyzed the for- 
mation history of two Milky Way like disk galaxies in 
detail and find gas accretion rates, and hence bolometric 
X-ray luminosities of the haloes, 6-7 times larger at z ~1 
than at z=0 for these disk galaxies. More generally, we 
find that gas infall rates onto these disks are nearly ex- 
ponentially declining with time, both for the total disk 
and the "solar cylinder" . This theoretical result hence 
supports the exponentially declining gas infall approx- 
imation often used in chemical evolution models. The 
infall time-scales deduced are ~5-6 Gyr, comparable to 
what is adopted in current chemical evolution models to 
solve the "G-dwarf problem" . 

It is moreover found for these two galaxies that galactic 
disks can form "inside-out" as well as "outside-in" , but 
in both cases the mean ages of the stars in the outskirts 
of the disks are >6-8 Gyr, broadly consistent with the 
findings of Ferguson & Johnson for the disk of M31. 

The amount of hot gas in disk galaxy haloes is con- 
sistent with observational upper limits. We "insert" the 
globular cluster M53 and the LMC in the haloes of the 
two Milky Way like disk galaxies and calculate disper- 
sion measures to these objects. Our results are consistent 
with upper limits from observed dispersion measures to 
pulsars in these systems. 

Our simulations indicate that we can reproduce the ob- 
served peak in the cosmic star formation rate at redshift 
z ~ 2. Depending on the star formation and feedback 
scenarios we predict either a monotonically decreasing 
cosmic star formation rate beyond these redshifts or a 
second peak at z^6-8 corresponding to the putative pop- 
ulation III and interestingly similar to recent estimates of 
the redshift at which the Universe was reionized. These 
various scenarios should hence be observationally con- 
strainable with upcoming instruments like JWST and 
ALMA. 

In this paper we have mainly been concerned with the 
present day properties of the simulated disk galaxies. We 
shall in forthcoming papers study the present day prop- 
erties of the E/SOs in detail as well as the dynamics of all 
types of galaxies. In particular, we intend to study the 
kinematics and dynamics of the stellar haloes of Milky 
Way like disk galaxies to compare with current models of 
the Galactic stellar halo (e.g., Sommer-Larsen et al. 1997; 
Helmi & White 1999; Helmi ct al. 1999). We will also 
study the higher redshift properties of all types of galax- 
ies and the properties of medium to high-z damped 
Lya systems predicted by our simulations. Moreover, 
we will incorporate chemical evolution including non- 
instantaneous gas and heavy element recycling in the 
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simulations and replace thermal energy by entropy as 
an independent variable in SPH (Springel & Hernquist 
2002). 

Last, but not least we will also study much more 
closely the physical factors determining the final Hubble 
type of the galaxies formed. It is, however, already 
at this point clear from our simulations that final 
galaxy morphologies primarily reflect the merging 
histories of the galaxies in the sense of the very de- 
tailed merging histories — just to give an example, 
we find that even major merging sometimes can be 
constructive in forming disk galaxies with extended 
disks. Another interesting result is that whereas the 
mean spin-parameter for the 5 spheroid dominated 
systems (E/SO) is < A 20 o >_e/ so =0.043±0.005, it is 
<A 2 oo>s=0.037±0.002 for the 35 disk galaxies. Hence 
our results do not seem to support the notion that the 
DM haloes in which disk galaxies form on average have 
larger spin-parameters than the haloes of early type 
galaxies (e.g., Mo, Mao & White 1998; van den Bosch 
1998). 
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